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ABSTRACT 

The simple, conventional dark matter halo mass definitions commonly used in cos- 
mological simulations ("virial" mass, FoF mass, M5o,ioo,20o,...) only capture part of 
the collapsed material and are therefore inconsistent with the halo mass concept used 
in analytical treatments of structure formation. Simulations have demonstrated that 
typical dark matter particle orbits extend out to about 90 per cent of their turnaround 
radius, which results in apocenter p assages outside of the cur rent "virial" radius on 
the first and also on the second orbit ( Diemand fc Kuhlen|2008 ) . Here we describe how 



the formation history of haloes can be used to identify those particles which took part 
in the halo collapse, but are missed by conventional group-finders because of their 
remote present location. These particles are added to the part of the halo already 
identified by FoF. The corrected masses of dark haloes are significantly higher (the 
median mass increase is 25 per cent) and there is a considerable shift of the halo mass 
function towards the Press & Schechter form. We conclude that meaningful quanti- 
tative comparisons between (semi-)analytic predictions of halo properties (e.g. mass 
functions, mass accretion rates, merger rates, spatial clustering, etc.) and simulation 
results will require using the same halo definition in both approaches. 

Key words: cosmology: dark matter - large-scale structure of Universe - methods: 
N-body simulations 



1 INTRODUCTION 

Theoretical models of cosmological structure formation are 
based on the ansatz that an overdense region in the linearly 
evolved density field collapse into a dark matter halo when 
the linear density contrast reaches a certain collapse thresh- 
old 5c- In an idealised, radial spherical collapse, a region 
collapses into a point at the time when its linear density 
contrast grows to 5c = 1.686 (|Gunn & Gott||1972|). |Press 



|fc Schechter] ( |1974[ PS hereafter) have used this collapse 
threshold to calculate the abundance of haloes from the mass 
variance o"^(M) of the linear density field. PS was extended 
by [Bond et aLl|T99]] and many others using the 'excursion 
set' approach, which allows a more robust derivation of the 
original PS mass function and also the prediction of further 
halo properties like accretion rates, merger rates and spatial 



clustering (see Zentner||2007 for a recent review of this ap- 
proach). Note that in the PS formalism the mass enclosed 
in the linear over-density and the resulting halo mass are 
assumed to be equal, i.e. the entire collapsing region makes 
up the final mass of the halo. 

Since the collapse itself is fairly complicated (i.e. non- 
linear, clumpy, non-radial, non-spherical) event, theoretical 
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models of structure formation cannot predict the detailed 
properties of the haloes they describe, i.e. PS does not de- 
scribe (or depend on) the mass distribution in their haloes. 
Modern cosmological simulations on the other hand calcu- 
late the z = non-linear density field reliably, but to ex- 
tract halo properties from simulations requires an opera- 
tional halo definition, i.e. the choice and implementation of 
some kind of halo finder. Unfortunately the conventional 
halo definitions used by simulators are still based on an out- 
dated, oversimplified picture: By imposing that a homoge- 
nous collapsing sphere reaches virial equilibrium, one finds 
that the material would settle within a sphere of half its 
turnaround radius, defined to be the virial radius rvir- By 
definition rvir contains the virial mass Mvir, which is as- 
sumed to be equal to the mass of the collapsed homogenous 
sphere. The moment of this virialisation is assumed to be 
when a radially collapsing sphere falls into one point. At 
this time the virial radius encloses Am = 178 times the 
mean matter density in a flat universe with Q.m = 1- This 
simple picture motivates the operational definition of haloes 
as spherical overdensities (SO hereafter) with p = 178pcrit 
|1992| . A related halo definition is the 



(e.g. Warren et al. 



friends-of-friends (FoF) algorithm (e.g. Davis et al. 1985 1 



FoF recursively links particles closer than some fraction h of 
the mean particle spacing. The conventional value of 6 = 0.2 
produces groups with comparable mean density as the SO 
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virial mass definition. The spherical collapse in ACDM leads 
to a larger Am and therefore a smaller linking length b (see 
White"2001', '2002; Wa rren et aLl[2006l [Tinker et aLl[2508l 
Robertson et al. 2009; for detailed comparisons of the halo 
mass function obtained with SO(Ai\/) and FoF(6) for a wide 
ranges in Am and b). 

SO and FoF are well defined and easy to use for the anal- 
ysis of large data sets from cosmological N-body simulations. 
SO halo masses correlate rather tightly with some galaxy 
cluster observables (SZ, X-ray, optical), especially for large 
Am values (see Tinker et al.||2008 and references therein). 
While FoF and especially SO based halo mass function are 
useful for comparisons with observations, they are not suited 
for comparisons with predictions from PS and related mod- 
els, because they only capture a fraction of the collapsed 
mass: Simulations show that spheres which enclose the final 
Mvir of a halo, collapse only by about a factor of 1.4, i.e. 
much less than the factor of 2 assumed in the definition of 
the virial radius ( Diemand, Kuhlen &: Madau|2007 1 . Related 
indications that haloes are much more extended than their 



2 THE CODE 



are their larger virialized regions (Maccio, Murante & 
Bonometto|2"003 Prada et al.|2006"l [Cuesta et al.|2008| ) aiid 



the fact that many haloes found between rvir and 2rvir of a 
large host halo are orbiting through and around this host, i.e. 



many were well inside its rvir at some earlier time ( 


Moore, 


Diemand & Stadcr'2004| |Gill, Knebe & Gibson 2005| [Lud- 


low et al.|2009). 


Diemand & Kuhlen 


( 2008p have shown that 



the dynamics in outskirts of haloes are very well described 
by the self-similar secondary infall model from [Fillmore fc| 



Goldreich (19841 and Bertschinger (1985): Typical particle 



(and subhalo) orbits extend out to about 90 per cent of 
their turnaround radius and the apocenter distance decays 
only very slowly and asymptotes to about 83 per cent of the 
turnaround radius. This implies that after falling into a halo 
a typical particle has its first two apocenter passages beyond 
the current rvir of this halo ( [Diemand fc Kuhlen|2008| . This 
causes significant amounts of is material which forms part 
of the collapsed region, and therefore of the halo mass in 
PS type models, to lie beyond rvir- The mismatch between 
the PS mass concept and Mvir is obvious and implies that 
quantitative comparisons between models and simulations 
will require that both approaches agree on what exactly is 
meant by 'halo'. 

One way out would be to combine analytic models with 
secondary infall to convert the total collapsed mass into a 
SO(Am) mass for comparison with simulation results. A sec- 
ond solution is to maintain the PS halo mass concept and 
to develop a group-finder which is able to extract the entire 
collapsed halo mass from cosmological simulations. We have 
investigated this second route and present in this paper an 
algorithm which finds the total mass of haloes in the PS 
sense, based on a standard groupfinder (FoF is used in this 
work, although SO provides similar results) and on the for- 
mation histories of haloes given by their merger trees. 

This paper is structured as follows. After this overview 
of the mass problem, we describe in Section [2] first the simu- 
lation setup we used for this work and thereafter the time de- 
pendent group finding algorithm. In Section [S] the resulting 
mass function based on this correction is presented whereas 
a summary and discussion is given in Section [4] 



In this section we present our correction code which com- 
putes the total halo mass. Sorting as well as searching algo- 
rithms are taken from Numerical Recipes ( Press et al.|1992 1 



2.1 N-body Simulations 

For this work we make use of two different ACDM simu- 
lations, whichwere run using the parallel tree code PKD- 



GRAV ( |Stadel|2001[ ). Both runs have a total number of 512'' 
particles, simulated in cubic volumes of 80/i~^ Mpc (Run 1) 
and 250/i~^ Mpc (Run 2) respectively, leading to particle 
masses of 3.18 x lO^h'^ Mq and 9.67 x lO^h'^ Mq. This 
setup enables us to probe a mass range of almost five orders 
of magnitude, although the statistics at the high mass end 
is not sufficient to give a detailed prediction for the cluster 
mass function. We are using a WMAP 1°' year cosmology 
for both simulations. The parameters are 



(Qq, ridm, ^A, h, (jg, n) 
(0.3,0.255,0.7,0.7,0.9,1) 



(1) 



The initial conditions for both simulations were created with 
GRAFIC2 ( Bertschinger|200T \ ; the analysis described in the 
following is based on 40 snapshots in time, ending at the 
present epoch. 



2.2 Finding the Total Halo Mass 

According to the PS ansatz, the mass of a dark matter halo 
is defined to be the entire mass which once collapsed from a 
connected (spherical or ellipsoidal) and overdense region. In 
a first attempt to find this mass, all particles in haloes today 
were traced back to their initial positions. Then we tried to 
find a well defined way to complete the initial particle dis- 
tribution into a convex and connected form. Unfortunately, 
these primordial regions turned out to be highly irregular 
and they are often not connected, so it was not possible 
to come up with a method to identify all particles which 
will collapse into a halo using only the initial distribution 
of those particles that ended up in haloes identified by the 
FoF groupfinder. We therefore consider the entire formation 
history of each halo in order to compute the total collapsed 
mass. The goal is to detect all particles which did belong 
to a progenitor of a redshift z = FoF-groujI^ and to join 
them together with the part of the halo already found by 
FoF into one complete z — halo. 

First, a halo merger tree is constructed based on the 
FoF group snapshots with a linking length b — 0.2. This 
is done by comparing the number of particles which went 
from a progenitor group G^i to a group Go at present time 
(A/gj.-jGoi the subscript denotes redshift z = 0), with the 
number of particles which were in the progenitor Gzi , 

f^G^.^G, ^2) 



and demanding that TZ ^ 0.5 Wzi to ensure unique remnants 
and merger tree without splits. 



^ We focus on 2 = 0. The same method can be applied at higher 
redshifts, if enough earlier snapshots are stored and analysed. 
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Figure 1. Projection on tiie x-z-plane of a cluster-sized tialo at 
z = 0. The left panel shows all particles which are members of 
the clusters z = FoF group (red), and the particles in a smaller, 
nearby FoF group (blue) . After correction for lost particles (right 
panel), many more particles, including the small FoF group now 
form part of the total mass of the cluster. The trajectory of the 
small group reveals that it is actually a subhalo (its irregular 
shape is due to tidal stripping, some of its tidal stream is captured 
by its FoF group). 



Now the total collapsed mass for each object is found 
by going through all particles and assigning each one to the 
correct z = halo. If a particle was part of a progenitor, 
the code assigns it to the unique successor according to the 
merger tree. Particles which do not already form part of the 
correct successor FoF group are added to that halo. Most of 
the corrected particles do not form part of any z = FoF 
group, i.e. they have orbited through a progenitor halo and 
are now found somewhere in the low density outskirts of the 
z = successor FoF group. Some of the corrected particles 
are members of a different 2 = FoF group, i.e. they form 
part of a subhalo which did fly through a progenitor halo 
and is now found outside the extent of the z = successor 
FoF group. In a conventional analysis such a subhalo would 
be misclassified and counted as an individual field halo. Our 
approach realises that this halo used to be a subhalo and 
should still be counted as a satellite, its mass is added to 
the primary halo and the subhalo disappears from the halo 
mass functiorj^ Note that such extended, but bound orbits 
are typical for particles and subhaloes which are on their 
first two orbits after infall ( [Diemand Sz Kuhlen||2008 1 . 

An illustration of one example outcome of our correc- 
tion method is shown in Fig. [T] In the FoF result (left 
panel) there are two distinct haloes, a cluster-sized halo (red, 
AfpoP = 1.16 X 10^'* Mq) and a nearby smaller halo (blue, 
M ~ lO" M0). Following the orbit of the small halo back 
in time shows that this halo is orbiting around the cluster 
and that it did form part of the clusters progenitors at ear- 
lier epochs, i.e. it is actually a subhalo which happens to 
lie outside its z = primary. Therefore our correction also 
takes care of all these subhaloes and adds them to their hosts 
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Figure 2. Phase-space distribution of dark matter particles at 
rcdshift ^ = of a small galaxy-sized halo. Particles in the FoF 
group are plotted as blue circles, red triangles give all halo mem- 
bers after the correction. The black dots are particles which are 
neither in the FoF group, nor in the corrected halo. In addition 
to the virial radius (solid line), the static radius (dot-dashed) is 
given ( [Cuesta et al.|2008| . Depending on the definition, the mass 
of this halo is Mvir = 2.13 X lO" Mq, Mp^p = 2.37 X lO" Mq, 
M.tatic = 4.73 X 10" M0 or Mtot^i = 3.26 x lO" Mq. 



in addition to the many other particles which were lost and 
not captured in the hosts z = FoF group. 

Fig. [2] shows the phase-space distribution of an FoF 
group and of the corresponding complete, corrected halo. 
The bulk (almost 90 per cent) of the additional particles 
from the correction are located between one and two virial 
radii. However, about 5 per cent of the corrected particles 
have travelled out to a distance of more than three virial 
radii after their infall. These exceptional distances suggest 
that these particles have gained kinetic energy during their 
passage through the halo, due to real (e.g. three body en- 
counters involving a massive subhalo. Sales et al.||2007 l or 
numerical effects. Some of these very remote particles might 
never fall back into this halo. At z = they might be un- 
bound or even bound to another halo. Nevertheless, for the 
total halo mass and especially for the purpose of compari- 
son with analytic results, we advocate counting even those 
remote particles as halo members, since they did fall in once 
and PS and related approaches do not account for the com- 
plicated nonlinear effects which led to their ejection]^ 



^ A particle can only be part of one halo, i.e. in our corrected 
mass function no mass element can be counted more than once. 



^ In a few cases two nearby 2 = FoF groups, which according to 
our merger tree do not appear to have ever been one and the same 
halo (i.e. merged and de-merged), might still have exchanged a 
significant number of particles. If the correction assigns more than 
15 per cent of the members of the smaller group to its larger 
neighbour we assume that the small group is a subhalo and add 
it to its primary halo. Such haloes did have a close encounter with 
the main halo, but they were not caught inside the main halo at 
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between Mstatic and Mtotai is that Mstatic includes all par- 
ticles within an increased search radius around each halo: 
If two FoF groups lie close together they do both prevail 
and grow in mass, and some of the additional mass may be 
assigned to both haloes. In the case of Mtotai every particle 
can contribute to at most one halo. If two nearby FoF groups 
emerge from a common progenitor, they are joined together 
into one larger halo. For galaxy-sized haloes our total-to- FoF 
mass ratios are smaller than the typical static-to-virial mass 
ratio in Cuesta et al. ( |2008| ), whereas it is the other way 
around for clusters. Our correction is slowly decreasing with 
mass (see linear fit in Fig.[3|, while the static-to-virial mass 
ratio is an increasing function at the low mass end. Above 
the characteristic mass (Af*) the typical static-to-virial mass 
ratio decreases quickly because the relative importance of in- 
fall increases, which pushes the static radius back close to 
the virial radius. Our method shows that haloes well above 
Af*, e.g. cluster haloes at z = 0, still do have a total mass 
which is significantly larger than their FoF and their virial 
mass. Below Ah mass corrections (up to 2-3 times the FoF 
mass) are more common, causing the 84* percentile to in- 
crease. 



Figure 3. The total collapsed mass-to-FoF mass relation at 2: = 
0. The solid line shows the median, the dot-dashed lines the 16"" 
and 84*'' percentile. The red dot-dashed line is a linear fit to 
the median: Mtotal/AfpoF = 1-25 - 9 X IQ-^ log(MpoF/3.16 X 
IQIOMq). 

Fig. |3] shows the total collapsed mass-to-FoF mass rela- 
tion for all dark matter haloes at redshift z = 0. For clarity, 
only the median as well as the IG*'^ and 84"^ percentiles 
are plotted. At the low mass end the scatter is largest: some 
haloes are up to three times heavier than their corresponding 
FoF group, while other haloes have about the same mass be- 
fore and after the correction. The median correction values 
are almost independent of mass, ~25 per cent for galaxy- 
sized haloes and ~20 per cent for clusters. The mass de- 
pendence is weak, because the secondary infall pattern with 
subsequent apocenter passages beyond the virial radius ap- 



plies to both galaxy- and cluster-sized haloes (Diemand & 



Kuhlen 20081. Only the weighting of these patterns differs 



somewhat: in clusters a larger mass fraction is currently in- 
falling and the relative importance of the older, distant ma- 
terial near its first and second apocenter is slightly smaller. 

The static mass (Cucsta ct al. 2008) is the mass of all 
particles enclosed in the static radius. By definition, the 
static radius is the largest radius where the mean radial ve- 
locity averaged over all haloes in this mass bin is still close 
to zero, i.e. < 0.05vcirc{rvir)- For the galaxy-sized halo in 
Fig. 2, the static mass is about 40 per cent larger than the 
our total mass. This originates from the fact that Afstatic en- 
closes all particles up to rstatic and many of those have never 
experienced infall (black dots in Fig. [2|. Another difference 



the analysed snapshots. An overlap of less than 15 per cent is 
assumed to be due to ejected particles and the two haloes are 
kept in the catalogue as individual objects. The threshold value 
of 15 per cent is arbitrary. However, we have checked that our 
resulting mass function is practically unaffected by this threshold 
value, as long as it is larger than a few percent. 



3 THE MASS FUNCTION 

3.1 Correcting the Mass Function 

In this Section we present the halo mass function (MF here- 
after) based on the total halo mass, instead of the conven- 
tional choice of FoF or SO halo masses. The PS ( [Press fc| 
Schechter|1974 1 result, as outlined in Section [ij is given by 

n{M,t)dM = ±-fps4^dM, (3) 



dlogA/ 



where p is the mean background density of the universe, 
" and 



fpsicr) = 



!^ CXp 



(4) 



is the multiplicity function which gives the fraction of mass 
associated with haloes in a unit range of logiy. |ShethfcTor-| 
( 1999 ST hereafter) proposed a more general formula 



motivated by the assumption of an ellipsoidal collapse, 
1 



/st(o-) = A 



"^11 + 

TV L 



u exp 



(5) 



where A = 0.322, a = 0.707 and p = 0.3 (see Desjacques 
|2008[ [Robertson et ar]|2009| for a discussion of elliptical 
collapse and collapse barriers). In the limit of A = |, a = 1 
and p = it reduces to the spherical collapse and to the PS 



form. Recent work by Bhattacharya et al. ( 2010) suggests a 
new fitting function, similar to the ST form, but with one 
additional parameter. Using 67 high resolution simulations, 
they found that ST deviates up to 40 per cent from their 
simulated FoF mass functions at the high mass end (Fig. 4 

their MF 



takes the form 



Bhattacharya et al. (20101). At redshift 



rST , s 



0.333 \/ — exp 

TV 



1 + 



0.788(5g 
2a2 



(6) 



0.78852 



0-807^ /<5c^/a788\ 1-795 
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Figure 4. Halo mass functions extracted at z = form our two 
simulation boxes, using the conventional FoF (b = 0.2) halo mass 
definition (green triangles and squares corresponding to Run 1 
and Run 2 resp.) and using the corrected total halo masses (red 
stars and circles corresponding to Run 1 and Run 2 resp.). The 
errorbars represent the Poissonian error. For comparison the PS, 
ST and [Bhattacharya et al.| | |2010[ l mass functions (black solid, 
blue dashed and cyan dotted lines) are plotted. The purple, dot- 
dashed line shows a fit (see Eq. [7| to our corrected mass function. 



The FoF mass functions from our two simulations, with 
a linking length b = 0.2, is plotted with green data points in 
Fig. |4] Only haloes with more than 50 particles are consid- 
ered, yielding a lower mass hmit of 2.3 x 10^° Mq. At the 
low mass end, the ST form of Eq. |5] provides a good fit, but 
it significantly starts to deviate at high masses. The recently 
proposed modified ST fitting formula of Eq. [6] fits our FoF 
based mass function very well over the entire mass range 
probed by our two simulations. 



[Bhattacharya et al.|(|2010[) as well as other empirical fit- 



ting functions (e.g. 


Jenkins et al.||2001 Warren et al.||2006 


Tinker et al. 2008 


Reed et al.||2007l are based on conven- 



tional mass definitions like Mvir, Mfop, M50, 100,200,..., which 
capture only a fraction of the total collapsed mass. They are 
therefore not suited for comparisons with analytical models 
of the PS and ST type. The corrected mass function (CMF 
hereafter), using the total halo masses after our correction, is 
plotted as red stars and circles in Fig. |4] Here the errorbars 
represent the Poissonian error. Below Af^ there is a shift of 
the mass function of about 15 per cent relative to the one 
based on FoF halo masses. In the cluster regime this shift 
seems to increase, but in order to get more reliable predic- 
tions a larger box simulation is needed. 

Based on the minimal-parameter multiplicity func- 
tion described in [Warren et al.| ( [2006[ ) and by using the 
Levenberg-Marquardt algorithm, we determine a best fit for 
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Figure 5. Residuals between the measured FoF(fe = 0.2) mass 
function and the CMF (solid line) at 2 = 0. In addition also 
the residuals between the measured mass functions analysed with 
varying linking length parameter, b = 0.27, 0.30, 0.33 and 0.40 
(dotted, dashed, long-dashed and dot-dashed), are shown. The 
left and right data sets correspond to Run 1 and Run 2 respec- 
tively. 

the CMF of the form: 

/cmf(o-, z = 0) = a{a'' + c) exp ^ - , 

a = 2.825, 6 = 0.138, c = -0.883, d = 0.154. 

(7) 

This fitting formula is represented by the purple line in Fig. 
II 

3.2 Comparison with Higher Linking Lengths 

As shown in Fig. |3| the total collapsed mass for individual 
haloes lies approximately 20 to 25 per cent above the the 
FoF(b = 0.2) mass. This is also refiected in the representa- 
tion of the mass function in Fig. [4] A simple way to increase 
the mass of a halo is to increase the linking length parameter 
in the FoF algorithm (or similarly relax the virial overden- 
sity criterion in an SO group finder). In this subsection we 
show that the total collapsed mass and therefore the CMF 
cannot be reconstructed by simply varying the only free pa- 
rameter in the FoF group finder. Fig. [5] shows the effects of 
increasing the FoF linking length and compares them with 
our CMF. None of the FoF mass functions with increased 
linking length are able to reproduce the CMF, even though 
the deviations between the CMF and FoF(6 = 0.27 - 0.33) 
mass function tend to get smaller above M*. 

For an even more unrealistic value of the linking length, 
b = 0.40, the deviations get quite large. In this case, other 
unphysical characteristics like interconnections of groups 
through slightly overdense filaments start to dominate and 
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Figure 6. Deviations from universality of the measured FoF 
(dashed lines) and CMF mass function (solid lines) at redshifts 
z = 1 (blue), 2 = 2 (red) and z = 3 (green). The horizontal lines 
represent the 20 per cent threshold often found in the literature. 
Correcting the halo mass as described in Section |2.2| shifts the 
mass function towards universality. 



therefore changing the overall shape of the mass function 
due a decrease of the the total number of galaxy-sized ob- 
jects and the sudden occurrence of super-clusters of order 



3.3 Universality of the Corrected Mass Function 

The extended PS formalism is from its analytic arguments 
expected to be universal, i.e. its predicted halo abundance 
depends only on the variance of the density field a{M), but 
not on redshift or cosmology (e.g. Zentner||2007 K However, 
it has been shown by various authors that universality only 
holds at the 10 to 20 per cent level when using FoF mass 
functions and slightly worse if haloes are identified with an 
SO group finder (e.g. |Lukic et al.|r2007| [Tinker et al.||2008 



Bhattacharya et al.|20l6 l 

To check if the CMF follows the universal behaviour, 
we ran two additional simulations with the same cosmology 
as in Eq. [l] but with a significantly higher initial redshift, 
and box sizes of 40/i~^Mpc and 80/i~'^Mpc respectively. By 
keeping the number of particles constant, these simulations 
provide a sufficiently high resolution at redshifts z = 1, 2, 3 
not only to recover the FoF mass function, but also to ap- 
ply our correction algorithm (i.e. tracing back particles over 
several time steps). For testing universality it is preferable 
to represent the mass function in the /(cr) — Incr"^ plane, 



M dn 



(8) 



po dlncr"^ ' 

since cosmology and redshift dependences are absorbed in 
(j(M). 



In Fig. [6] the deviations from universality of the mea- 
sured FoF and CMF mass functions at redshifts z — 1,2 and 
3 are shown. Whereas the FoF mass function is universal 
approximately at the 15 per cent level up to redshift z — 3 



(consistent with previous studies (e.g. Bhattacharya et al 



2010 and references therein)), we notice a significant shift 



towards universality in the CMF case for redshifts z > 1. 
This should not be surprising since the goal of the CMF 
was to approach the halo mass definition in the sense of PS 
as close as possible. 



Courtin et al. ( 2011 1 investigated the behaviour of mass 



functions in varying dark energy models. They found that 
the virialization process contributes to shaping the mass 
function in a cosmological and redshift dependent way and 
therefore a clear break down of universality when dark en- 
ergy starts to dominate {z < 1). In the case of the total 
collapsed mass, the virialization process is made more inde- 
pendent of redshift and hence the CMF is expected to show 
a higher degree of universality. 



4 SUMMARY AND DISCUSSION 

In this letter we present a group finding algorithm which, 
based on a conventional FoF group finder, identifies the to- 
tal collapsed mass of every halo in the sense of the PS the- 
ory and its extensions. We found a considerable shift of the 
mass function of order 15 per cent towards higher halo abun- 
dances, i.e. towards the PS form. The median mass increase 
for individual dark matter haloes is of order 25 per cent for 
galaxy-sized objects and of order 20 per cent for clusters. 
However there is still a considerable gap between our cor- 
rected mass function (Eq. [7| and the PS mass function (Fig. 
|4|. This can be partly understood by again tracing back all 
particles, that form the total mass, to their initial positions. 



Many of the disconnected FoF regions (see Section 2.21 are 
now connected, but not all of them. The typical initial region 
which will form part of the total mass of a 2: = halo still 
differs form a simple spherical or ellipsoidal overdensity, even 
though the difference does decrease when the total mass is 
used instead of the FoF or SO mass. Furthermore since the 
total mass correction is more or less uniform over the en- 
tire mass range, we find a higher abundance of clusters and 
superclusters than predicted by PS. This suggests that the 
PS ansatz might be too simple to allow precise predictions 
about abundance and properties of dark matter haloes. 



The classical secondary infall model | Fillmore & Goldreich 



1984 Bertschingcr 1985) describes the median orbits of par- 



ticles and subhalos in the outer regions of halos very accu- 
rately ( Diemand fc Kuhlen||2008 1. Extensions of this model 
match the mass distributions in and around simulated halos 
( [Ascasibar, Hoffman &c Gottlober||2007[ ). It could therefore 
be feasible to calculate the expected halo mass fraction be- 
yond the virial radius within this framework and it would 
be useful to compare such results with our measured median 
corrections (Fig. 3). 



Robertson et al. ( 2009 1 used cosmological simulations 



to test the excursion set ansatz by identifying the locations 
in the linear overdensity field that later collapse to form dark 
matter haloes. They found an inconsistency between the ef- 
fective collapse barrier of simulated haloes and excursion 
set formalism predictions for their abundance, and conclude 
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that the excursion set ansatz fails, i.e. that the extended PS 
formaUsm cannot predict halo abundances exactly. However 
only the common mass definitions (FoF and SO) were con- 
sidered and it would be worthwhile to reinvestigate these is- 
sues using the total collapsed mass instead. Using the total 
mass will also allow to derive physically meaningful mass 
accretion and merger rates from cosmological simulations, 
and it will change halo formation times and their spatial 
clustering. 

Taking the total mass of a halo into account, instead 
of adopting an artificial truncation at the virial radius, does 
make a difference in the interpretation of the kinematics of 
groups and clusters. For example, it increases escape veloci- 



ties significantly. The fast neighbours AndXIV ( |Majewski et 
aL]|2007| ) and AndXII ( [Chapman et al.|2007| ) might well be 
bound to Andromeda without requiring a much larger virial 
mass. And even objects which are quite isolated today, e.g. 
Tucana ( |Mateo|[T998| , might well have had a close interac- 
tion with their primary (in this case the Milky Way) at some 
earlier time and might have another one in the future. 
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